PrecSezione 9.7: Unità di elaborazione dei segnali SuCapitolo 9: Densità spettrale e filtraggio Sezione 9.9: Filtri analogiciSegue

## 9.8Filtri digitali↓

Sono una particolare classe di filtri, per i quali l’integrale di convoluzione si riduce in una sommatoria. Sebbene nel seguito sia analizzato il solo caso di segnali tempo-continui, i filtri digitali sono di particolare rilevo perché permettono di svolgere le operazioni di filtraggio operando direttamente sui campioni dei segnali (vedi § 4.5.1↑), e dunque possono essere realizzati via software o hardware logico [445] [445]Il tema delle realizzazioni numerihe dei filtri digitali non è al momento qui svolto, ma citiamo come fonte di approfondimento `http://www.dspguide.com/ch14/6.htm`)..

### 9.8.1Filtro trasversale ↓

E’ una particolare architettura di filtro digitale, mostrata nella parte sinistra di fig. 9.16↓, e per la quale è facile verificare come la risposta impulsiva abbia espressione
(17.34) h(t)=N⎲⎳n= 0cnδ(t−nT)
in cui N è l’ordine del filtro [446] [446]I coefficienti cn vengono indicati nei testi inglesi come taps (rubinetti) in quanto possono essere pensati “spillare” frazioni del segnale. Per effetto di un processo di trasposizione linguistica, gli stessi coefficienti in italiano vengono a volte indicati discorsivamente come tappi (!)., che essendo finito, fa rientrare il filtro trasversale nella categoria dei filtri fir↓ (finite impulse response). E’ altrettanto facile verificare che la risposta in frequenza associata alla (17.34↑) risulta
H(f)= ℱ{h(t)}=N⎲⎳n= 0cne−j2πfnT
Una tale architettura può essere rappresentativa di un effettivo fenomeno naturale, come ad esempio la presenza di “echi” nel segnale ricevuto. Come mostrato a destra in fig. 9.16↓, H(f) risulta periodica (in frequenza) con periodo F=(1)/(T), dato che tutti gli esponenziali e−j2πfnT lo sono. Pertanto questa architettura può essere adottata per sintetizzare una risposta in frequenza desiderata in una certa gamma di frequenze, e “qualsiasi” altrove. Infatti, se il segnale di ingresso x(t) è limitato in una banda (base o traslata) di estensione minore di (1)/(T), l’azione filtrante ha luogo appunto nella sola banda del segnale.
Figura 9.16 Schema simbolico di un filtro trasversale e sua risposta in frequenza
Per sintetizzare il filtro a partire dall’andamento desiderato di H(f) nella banda di interesse, si calcolano i coefficienti ck mediante la formula
(17.35) ck=T1 ⁄ 2T⌠⌡− 1 ⁄ 2TH(f)ej2πfkTdf
che è del tutto analoga alla eq. (6.21↑) di pag. 1↑ od (a parte il segno di f) alla eq. (6.2↑) di pag. 1↑. Ovviamente, se H(f) è qualsiasi, occorrerebbero un numero infinito di coefficienti ck; usandone un numero inferiore (finito) si produce una approssimazione
della H(f) desiderata [447] [447]Il troncamento (§ 2.2.2↑) della serie di coefficienti ck avviene in modo simmetrico rispetto a c0, prendendo cioè sia gli indici positivi che quelli negativi. Viceversa, nello schema di filtro trasversale si usano solo coefficienti con indici ≥ 0. Nel caso in cui l’H(f) da cui partiamo sia reale pari, allora i ck sono una serie reale (pari), garantendo un filtro idealmente realizzabile, ma la cui h(t)=∑N⁄ 2k=  −N⁄ 2ckδ(t−kT) necessita di una traslazione temporale per essere anche fisicamente realizzabile. Se invece H(f) ha un andamento qualunque, non si può dire nulla a riguardo di eventuali simmetrie per i coefficienti ck..

### 9.8.2Filtro trasversale del 1o ordine

E’ descritto dalla architettura mostrata a lato, a cui corrisponde una risposta impulsiva
h(t)=δ(t)+αδ(t−T)
la cui trasformata è H(f)= 1 +αe−j2πfT( [448] [448]In questo caso pur risultando H(f) a simmetria coniugata (H(f)=H*(−f)), è complessa. Pertanto, i coefficienti ck ottenibili dalla (17.35↑) sono reali, ma non necessariamente pari. Svolgendo i calcoli, si ha: ck=T∫1  ⁄ 2T− 1 ⁄ 2T(1 +αe−j2πfT)ej2πfkTdf=T∫1 ⁄ 2T− 1 ⁄ 2Tej2πfkTdf+αT∫1 ⁄ 2T− 1 ⁄ 2Tej2πf(k− 1)Tdf. Il primo integrale è nullo per k≠ 0, mentre il secondo per k≠ 1, in quanto le funzioni integrande hanno media nulla sull’intervallo 1 ⁄T; pertanto c0= 1 e c1=α, esattamente come è definita la risposta impulsiva.), e quindi
|H(f)|2 = (1  +αcos2πfT)2+(αsin2πfT)2= = 1 + 2αcos2πfT+α2(cos22πfT+ sin22πfT)= = 1 +α2+ 2αcos2πfT
A lato è riportato l’andamento di |H(f)|2 per due valori di α= ±.5, ed osserviamo che nell’intervallo di frequenze |f|<(1)/(2T) può assumere un comportamento passa-alto oppure passa-basso [449] [449]In altre parole, l’andamento ondulatorio di |H(f)|2 rende il filtro idoneo a diversi utilizzi, in funzione dell’andamento in frequenza del segnale di ingresso., in funzione del segno di α. Notiamo inoltre che ponendo α=  − 1 si ottiene un differenziatore, in grado di rimuovere dall’ingresso segnali periodici di periodo T.

### 9.8.3Filtro a risposta impulsiva infinita del 1o ordine↓

L’architettura rappresentata in figura identifica un filtro di tipo ricorsivo, ovvero il cui valore in uscita dipende dalle uscite precedenti. La risposta impulsiva corrispondente ha una durata infinita [450] [450]In questo caso si parla di filtro iir o infinite impulse response., ed è pari a
h(t)=∞⎲⎳n= 0αnδ(t−nT)
La funzione di trasferimento risulta quindi pari a H(f)=∑∞n= 0αne−j2πfnT e, ricordando la formula dello sviluppo in serie geometrica ∑∞n= 0βn=(1)/(1 −β), si può scrivere
H(f)=∞⎲⎳n= 0(αe−j2πfT)n=(1)/(1 −αe−j2πfT)
Notiamo subito che il filtro è stabile purchè |α|< 1, altrimenti si ha una uscita anche senza ingresso, ovvero uscita infinita con ingresso limitato. Per ciò che riguarda |H(f)|2, otteniamo
|H(f)|2=(1)/((1  −αcos2πfT)2+(αsin2πfT)2)=(1)/(1 +α2− 2αcos2πfT)
Le curve mostrate appresso rappresentano 10log10|H(f)|2, con T= 1 e diversi valori di α, positivi a sinistra e negativi a destra. Osserviamo infine che il caso α=  1 corrisponde ad avere un integratore perfetto che, ad esempio, produce una rampa in uscita, se in ingresso c’è un gradino.
PrecSezione 9.7: Unità di elaborazione dei segnali SuCapitolo 9: Densità spettrale e filtraggio Sezione 9.9: Filtri analogiciSegue